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Prograde rotation of protoplanets by accretion of pebbles 
in a gaseous environment 
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ABSTRACT 

We perform hydrodynamical simulations of the accretion of pebbles and rocks onto 
protoplanets of a few hundred kilometres in radius, including two-way drag force cou- 
pling between particles and the protoplanetary disc gas. Particle streams interacting 
with the gas far out within the Hill sphere of the protoplanet spiral into a prograde 
circumplanetary disc. Material is accreted onto the protoplanet due to stirring by the 
turbulent surroundings. We speculate that the trend for prograde rotation among the 
largest asteroids is primordial and that protoplanets accreted 10%-50% of their mass 
from pebbles and rocks during the gaseous solar nebula phase. Our model also offers a 
possible explanation for the narrow range of spin periods observed among the largest 
bodies in the asteroid and trans-Neptunian belts, and predicts that 1000 km-scale 
Kuiper belt objects that have not experienced giant impacts should preferentially spin 
in the prograde direction. 

Key words: hydrodynamics - Kuiper belt - minor planets, asteroids - planetary 
systems: protoplanetary disks - planets: rings - Solar system: formation 



1 INTRODUCTION 

Planets form in protoplanetary gas discs as dust grains 
collide and grow to ever larger bodies (ISafronovl Tl969l: 
iDominik etafl 120071 1 . Accumulation of differentiated me- 
teorite parent bodies must occur while 26 Al and 60 Fe arc 
still present in the solar nebula, constraining their forma- 
tion time to at most a million years after the formation 
of calcium-aluminium-rich inclu sions (|Bizzarro et al.l [20051 : 
lYang. Goldstein, fe Scottl 120071 '). This coincides with the 
gaseous solar nebula epoch and indicates that planetesimals 
and protoplanets of several hundred kilometres in size were 
accumulated in a dense gaseous environment. 

The alignment between the orbits of the planets and 
the rotation of the Sun is one of the key pieces of evidence 
for the nebular hypothesis of the formation of the solar sys- 
tem. The large majority of satellites also orbit their host 
planet in what is termed the prograde direction, i.e. orbit in 
the same direction as the planet orbits the sun, suggesting 
that they formed in circumplanetary discs. This regularity 
seems to extend to the predominantly prograde rotations of 
the planets and the largest asteroids. While the gas giant 
planets have likely been spun up by gas accretion, growth 
of solid planetary bodies by accretion of km-sized planetes- 
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imals is not generally expected to lead to systematic pro- 
grade rotation. Post-accretion, giant impacts which can sig- 
nificantly change the obliquity and spin rate of p rotoplanets 
(|Dones fe Tremainel Il993l ; iKokubo fe Idall2007l ) do so in a 
stochastic way yielding no preference for prograde motion. 
Successful attempts to explain a systematic prograde con- 
tribution from accretion of solids in a gas-free medium re - 
quire special surface density pr ofiles ( Lissauer fe Karv| | l99lf ) 
or "semicollisional" accretion (|Schlichting fe Sarill2007l ). i.e. 
collisions between small particles and subsequent accretion 
into a prograde circumplanetary particle disc that feeds the 
central object. 

In this paper we show that prograde rotation may be 
a natural consequence of planetesimal growth from accre- 
tion of cm-sized solids in a gaseous environment. Providing 
a physical mechanism for prograde rotation in turn supports 
that the dominance of prograde rotation among the largest 
asteroids is significant despite their low number. Our model 
assumes simultaneous presence of cm-sized "pebbles" and 
gas. Such conditions are supported by observations of pro- 
toplanetary dis cs, which show evidence of dust growth to a t 
least a few cm (|Testi et al.ll2003l ; ISicilia-Aguilar et al.ll2007l V 
The inferred population of mm-cm pebbles ca n be very sub- 
stant ial, on the order of 10~ 3 solar masses l|Wilner et al.l 
l2005h . 

We present hydrodynamical computer simulations of 
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Figure 1. Accretion of pebbles onto a protoplanet fixed in the centre of the coordinate system. The first panel indicates the Hill radius 
and the accretion radius (inner hole). The colour coding shows the local density of solids, in units of the mid-plane gas density. Particles 
are initially accreted from the second and fourth quadrants (second panel). As pebble streams drag the gas down towards the protoplanet, 
the incompressible gas expands in great particle-free bubbles, eventually making the system turbulent (third panel). A prograde accretion 
disc, fed from the second and fourth quadrants, forms around the protoplanet (fourth panel). 



the accretion of pebbles and rocks onto young protoplan- 
ets in a gaseous nebula environment. We find that parti- 
cle streams interacting with the gas far out within the Hill 
sphere dissipate enough energy that the particles spiral down 
towards the protoplanet and form a prograde accretion disc. 
The prograde motion is transferred to protoplanet rotation 
as the particles are accreted. To explain the observed rota- 
tions of the largest asteroids and Kuiper belt objects, we pro- 
pose that protoplanets accreted a major fraction (10%-50%) 
of their mass from pebbles and rocks during the gaseous so- 
lar nebula phase. 

The paper is organised as follows. In fj2]we describe the 
set up of our simulations. In f|3] the main results are pre- 
sented, both for 2-D and 3-D simulations. We speculate in 
fj4] about implications for the observed prograde rotation of 
large asteroids and the future prospect of measuring rota- 
tion directions of Kuiper belt objects. We conclude on our 
findings in Sj5] The appendices contain details that do not 
appear in the main paper and several tests of the validity of 
our results. In Appendix[X]the dynamical equations are pre- 
sented. We validate the code in Appendix [B] against known 
results on the accreted angular momentum. Appendix [C] de- 
scribes drag forces in more detail and also the translation 
from dimensionless friction time to physical particle size. 
In Appendix [D] we compare time-scales of drag force and 



particle collisions and conclude that drag forces are by far 
the dominant force, except in regions that are highly par- 
ticle dominated. In Appendix |E] we present results at lower 
and higher resolution and conclude that prograde rotation 
is strengthened with increasing resolution. The dependence 
of the results on the size of the simulation box are treated 
in Appendix [F] 



2 SIMULATION SET UP 

We perform numerical simulations of the coupled motion of 
solids and gas in a local coordinate frame corotating with the 
disc at the Keplerian freq uency Ok at an arbitrary dis tance 
ro from the central star (|Goldreich fc Tremaineiri980l ). Co- 
ordinate axes are oriented such that x points radially away 
from the star and y points along the rotation direction of the 
disc. Particles are treated by a Lagrangian method, interact- 
ing with the gas through two-way drag forces operating on 
the friction time-scale Tf. An additional gravitational force 
contribution from a protoplanet of mass M p fixed in the cen- 
tre of the coordinate frame is added. The gas is evolved on 
a fixed grid and also feels the gravity of the protoplanet (see 
Appendix [X] for details on the dynamical equations). In Ap- 
pendix [B] we validate the results of the code against known 
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Figure 2. Delivery of prograde angular momentum to a protoplanet of 400 km radius (if placed at ro = 10 AU). The dots show the 
angular momentum of a subset of the accreted particles, while the blue (black) and red (grey) lines indicate the instantaneous and 
cumulative average. (A-B) Accretion of pebbles from environment with mean solids-to-gas ratio e = 1.0 and e = 10.0, respectively. The 
presence of a Keplerian particle disc is clearly visible at late times where the accreted angular momentum is exactly equal to that of a 
Keplerian disc (dashed line). (C) In the low solids-to-gas ratio case (e = 0.1) the pebbles have little influence on the gas and their motion 
is strongly damped as they approach the protoplanet, imparting only little rotation to the protoplanet. (D) Without drag forces these 
effectively very large particles give a net retrograde rotation to the protoplanet. 



resul ts on angular momentum accretion l|Schlichting fc SarH 
120071 ). 

Momentum conservation in the two-way drag force cou- 
pling between solid particles and gas is assured by assigning 
the drag force added on a particle back to the gas with 
the opposite sign by a second order scheme. This drag force 
scheme has been shown to successfully reproduce the linear 
growth rate of streaming in stabilities in protoplanetary discs 
i|Youdin fc Johansedl2007l ). 

We integrate the dynamical equations of particles and 
gas using the high o rder finite difference code Pencil Code 
l)Brandenbura 12003). To avoid underresolution of motion 
near the protoplanet, we always freeze the gas inside a ring 
of diametre eight grid points around the protoplanet. Par- 
ticles crossing that ring are considered to be accreted onto 
the protoplanet. 

The Hill radius, R H = (GM p /3f2l) 1/3 where G is the 
gravity constant, bounds the region where the gravity of the 
protoplanet dominates over the tidal force from the star. 
In a protoplanetary disc with vertical pressure scale height 



H, a protoplanet with mass GM P = 10~ ti n^H 3 has 7? H « 
0.0071/. At r = 10 AU, with disc aspect ratio H/r = 0.06, 
the protoplanet has mass M p ~ 4 x 10 23 g and radius R p ~ 
400 km. At r = 3AU (disc aspect ratio H/r = 0.043) the 
protoplanet mass would be M p w 1.6 x 10 23 g and its radius 
R p w 290 km. 



We use a standard box size of (L x ,L y ) = 
(0.047/, 0.04//), at a standard resolution of 512 2 grid points 
and 10 6 particles. This gives a resolution element of around 
5,000 kilometres, making the inner cavity 20,000 kilome- 
tres in radius at the standard resolution. In Appendix [E] 
we present 256 2 and 1024 2 simulations. The inner cavity 
corresponds to 1/10, 1/20 and 1/40 of the Hill radius for 
low, medium and high resolution studies, respectively. The 
convergence test in Appendix [E] shows that the prograde 
rotation is strengthened at higher resolution with a corre- 
spondingly smaller inner accretion radius. 
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3 RESULTS 

In Figure[T]we show results of a simulation with an assumed 
solids-to-gas ratio e of unity (mimicking conditions in a sed- 
imentary mid-plane layer) and a friction time J?KTf = 0.1 
(relative to the inverse Keplerian orbital frequency). This 
friction time represent s particles of a few cm in radius 
(| Weidenschillind 1 1977al ) . See Appendix \C\ for details. The 
total mass of particles in the simulation domain is approxi- 
mately 1.8 x 10 22 g at r — 3 AU (with particle column den- 
sity S p = 3 gem" 2 ) and 6.4 x 10 22 g at r = 10 AU (with 
S p = 0.5 gem -2 ). 

As time progresses a spiral density wave structure ap- 
pears inside the Hill sphere. Particles are accreted primarily 
from the second and fourth quadrants. As the particles drag 
the gas along towards the protoplanet surface, and the gas 
is nearly incompressible, the gas leaves through a "wind" in 
the first and third quadrants. This wind prevents particles 
from entering from these sides. 

Later in the simulation the particle-gas layer becomes 
turbulent, due to the bubbles of gas expanding from the pro- 
toplanet surface. A prograde particle disc forms around the 
protoplanet, still fed from the second and fourth quadrants. 
This disc is completely dominated by particles (p p /p g 3> 10) 
and accretes slowly, due to stirring by the turbulent sur- 
roundings. We show in Appendix [D] that while particle col- 
lisions are negligible in most of the simulation domain, col- 
lisions will dominate over gas drag in the circumplanetary 
disc. However, we do not model the collisional viscosity in 
the presented simulations, as we wish to keep the system as 
simple as possible. A numerical algorithm for treating colli- 
sions between superparticles is currently under development 
(Johansen, Lithwick, & Youdin, in preparation). 

In Figure [2] we plot the angular momentum of accreted 
particles as a function of time, for a subset of randomly 
sampled particles. Initially, during the spiral density wave 
phase, there is a preference for accreting prograde particles. 
This is normally followed by a short phase of retrograde ac- 
cretion, where a retrograde accretion disc may temporarily 
form. Eventually the particles always enter a prograde ac- 
cretion disc, and the protoplanet accretes the full Keplerian 
surface frequency. 

Convergence tests are presented in Appendix [E] Both 
lower and higher resolution studies result in a similar 
amount of prograde accretion, although there is a trend for 
accreting stronger prograde rotation with increasing reso- 
lution. In Appendix [F] we present results for a box size of 
0.1H x 0.1H. The larger box size has a net positive effect 
on the prograde rotation of the protoplanet. The sign of the 
angular momentum is more consistently positive than in the 
small box, and there are fewer periods of retrogade accre- 
tion. 



3.1 Dependence on particle size 

It was shown analytically bv lMuto &: Inutsukal l|200ct ) that 
particles approaching the protoplanet feel the opposing ef- 
fects of protoplanet gravity (at large distances) and grav- 
itational scattering (at short distances). For particles with 
ftnTf < 1, drag forces are so strong that particles are ac- 
creted into the Hill sphere, while larger particles are scat- 
tered by the planet and only approach the protoplanet to 
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Figure 3. Accretion of larger and smaller particles. Marginally 
coupled particles i?KTf = 1 particles (top plot) readily form a 
prograde particle disc. The particles are much less affected by 
the gas than the particles presented in Fig. [2] thus the particle 
disc accretes very slowly. Smaller particles with f2K T f = 0.01, 
corresponding to sizes of a few mm, are so coupled to the gas that 
accretion takes place as a sedimentation towards the protoplanet, 
with a significant loss of angular momentum to the gas (bottom 
plot). 



some minimum distance. Similarly Wcidcnsc hilling fc Davis! 
|l985h showed that planetesimals above a certain size are 
trapped in mean motion resonances with growing planets. 
This shepherding might affect boulders of sizes metres to 
ten metres. 

In Fig.|3]we show the accretion of larger, J?Ki~f = 1, and 
smaller, f?KTf = 0.01, particles, corresponding to a few me- 
tres and a few millimetres, respectively. The larger particles 
readily form a prograde particle disc around the protoplanet. 
However, the marginal coupling to the gas makes the par- 
ticle disc much more stable to dynamical shaking by the 
surrounding gas, and the disc only accretes slowly. Smaller 
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particles, on the other hand, are very coupled to the gas, and 
accretion takes place as a slow sedimentation. The particles 
lose a significant fraction of their angular momentum to the 
gas and arrive almost perpendicular to the protoplanet sur- 
face. 

Thus we can put a relatively tight constraint on the 
particle sizes that can lead to prograde rotation of pro- 
toplanets. A dimensionless friction time between approxi- 
mately 0.1 and 1.0 corresponds to particles of radii between 
a few cm and a few ten cm. However, at several ten AU 
from the central star, where the gas column density drops to 
Eg m 5 gem -2 , the "good" particle sizes are between a few 
mm and a few cm, approximately corresponding to the peb- 
bles t hat are observed in the outer regions of protoplanetary 



bles that arc observed m the outer regions or protoplanetary 
discs jTesti et al.ll2003l; IWilner e t alj|2005l ; iRodmann et all 
2006 ; iLommen et al.ll2009I ). 



3.2 Sub-Keplerian motion 

An additional effect that we have not considered so far 
is the sub-Keplerian motion of the gas. As the gas is 
slightly pressure supported in the radial direction, the gas 
orbital motion is slower (by approximately 5% of the sound 
speed or 0.5% of the Keplerian speed) than the Keplerian 
value (jAdachi. Havashi. fc Nakazawal ll976; Wci denschillind 
Il977al ). Isolated particles face a constant headwind of the 
gas, causing them to orbit slower and to drift radially. The 
stationary protoplanet is thus exposed to a "sandblast" with 
velocity ~ 25ms -1 . In Fig. [4] we consider a simulation of 
i?KTf = 0.1 particles with an initial solids-to-gas ratio of 10. 
The strong mass loading of particles reduces the head wind 
significantly in the beginning, and we observe the usual pro- 
grade accretion. However, as the particles lose influence on 
the gas, a sandblasting phase occurs in which no net rota- 
tion is induced on the protoplanet. Overall the protoplanet 
nevertheless accretes a significant fraction of the Keplerian 
surface frequency. 



3.3 3-D results 

For simplicity, and to obtain higher resolution and smaller 
accretion radius, we have so far presented results of 2-D 
simulations. In this section we show that we obtain similar 
results in 3-D as in 2-D. We ignore the effect of vertical grav- 
ity and stratification on the gas, a small effect in boxes with 
a vertical extent of L z /2 = 0.02//. We initialize the particles 
with a Gaussian density distribution around the mid-plane, 
of scale height H p = 0.01H, and a column density of 0.01 
relative to the gas. The results of a 3-D simulation with 
256 3 grid points and 8 x 10 6 particles are shown in Fig. [5] 
The particles show the usual strong preference for prograde 
accretion. A prograde accretion disc also forms around the 
protoplanet in the 3-D simulation (see Fig.[6|, but this is not 
as evident from the accreted angular momentum as in 2-D 
because isolated particles, now entering from all directions, 
are accreted faster than the disc material. We conclude that 
our assumption of 2-D flow represents the true 3-D physical 
behaviour of the system well. The 2-D assumption, at the 
same time, allows for much higher resolution investigations 
and a higher particle number per grid cell. 
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Figure 4. Sandblasting a protoplanet with pebbles. (A) Mean 
particle density in the box as a function of time, (B) the accreted 
angular momentum. At the initial solids-to-gas ratio of 10, a net 
positive angular momentum is accreted by the protoplanet. But as 
the particles are accreted, the remaining particles drift stronger 
and stronger relative to the protoplanet, eventually leading to 
ballistic accretion with no net angular momentum. 
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Figure 5. 3-D accretion. The 3-D simulation has resolution N x = 
N y = N z = 256 and jV p = 8 X 10 6 . We recover the usual preference 
for prograde accretion, showing that the higher resolution 2-D 
simulations presented in the Fig. [2] represent the physical system 
well. 
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Figure 6. 3-D structure. The figure shows the column density of particles in the x-y plane (left panel) and in the x-z plane (right panel). 
The column density has been normalized by the gas column density in the relevant direction. A thin particle disc orbits close to the 
protoplanet, surrounded by the thicker particle mid-plane layer outside the Hill sphere. 



4 DISCUSSION 

4.1 Prograde rotation in the solar system 

The asteroid belt is a collisionally evolved collection of 
rocky and icy bodies residing roughly between the orbits 
of Mars and Jupiter. Analytical and numerical studies in- 
dicate that the largest asteroids have not suffered enough 
collisional evolution to s i gnificantly alter their spin vector 
l|McAdoo fc Burns! [l973l ; iFarinella et all ll992l V These ob- 
jects have thus likely retained their primordial spin angular 
momentum and may be re presentative of the orig inal proto- 
planet population (see also iMorbidelli et al.ll2009 ). Figure [7] 
depicts a collection of spin axis obliqu ities of 62 asteroids 
large r than 50 km in diametre (from iKrvszczvriska et al.l 
I2007I ). The obliquity is here defined as the angle between 
the pole vector of an asteroid and the normal to the eclip- 
tic plane. While smaller asteroids have collisionally random- 
ized obliquities, the largest ones (r > 150 km) display a 
propensity towards prograde rotation (obliquity < 90°), i.e. 
rotation in the same direction as the orbit. Indeed, the two 
most massive asteroids, Ceres and Vesta, spin progradely 
with obliquities of 9° and 36° and periods 9.1 and 5.3 hours. 
If their spin orientations were collisionally randomized, the 
probability that their mean obliquity is less than the ob- 
served value is only p ~ 2.5 x 10~ 3 (3<r). In Fig. [8] we show 
the mean pole angle (relative to the ecliptic) of asteroids as 
a function of their size. There is a trend towards prograde 
rotation for all size ranges, although the statistical signifi- 
cance is only 1-2 a. 

The Kuiper belt, a dynamically stable region beyond 
Neptune, p reserves roughly a dozen known 1000 km-scal e 
icy bodies (|jewittl [1993 : IChiang et alj|2007l : lBrownll200el . 
As in the case of asteroids, collisional evolution simulations 
of Kuiper belt objects (KBOs) show that bodies larger than 
r ~ 20 km have sustained re latively mild collisional evo- 
lution (|Davis fc Farineilalll99^ and their primordial spins 
have survived the ra ndom effect o f impacts for the age 
of the solar system (|Lacerda|[2005l '). No KBO obliquities 
have been derived as the known objects have moved less 
than 30° along their orbits since discovery, but several ro- 
tation periods have been measured (Fig. [9}. With the ex- 
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Figure 7. Obliquity of asteroids larger than 50 km in radius with 
respect to the ecliptic plane (x-axis). Each asteroid is represented 
by a blue (black) dot, connected to the origin by a grey line. The 
angle between each grey line and the positive y-axis is the obliq- 
uity of the corresponding asteroid. The length of the lines gives 
the spin frequency of the asteroid in rotations per day. Asteroids 
in the upper quadrant are prograde and asteroids in the lower 
quadrant are retrograde. Asteroids larger than r = 150 km, rep- 
resented by larger, red (dark grey) dots, show a preference for 
prograde rotation. 



ception of Pluto and Ha umea, both believed to have su f- 
fered massive collisions (|Canud 120051 ; iBrown et al.l 120071 ) , 
the rotation periods of KBOs large r than r = 200 km 
lie in a narrow 5 to 20 hour ran ge ([Lacerda fc Luull2006l : 
ISheppard. Lacerda. fc Orti3 120081 ) . The same is true for 
Ceres and the other three asteroids in the same size range, 
which have orbital timescales one to two orders of magni- 
tude shorter than KBOs. Such uniformity in spin period for 
bodies with different compositions and presumably formed 
and evolved at very different heliocentric distances suggests 
a common spin-up mechanism. 
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Figure 8. Mean pole angle of asteroids as a function of their size 
in kilometres. The black dots show the cumulative mean angle 
from the ecliptic. We connect the dots with a line to guide the eye. 
The first dot from the left is the pole angle of Ceres. The second 
dot is the mean of Ceres and Pallas, and so on. The shaded region 
is the 1-cr region expected from random spin orientations. The 
line is above the grey area until approximately the 10th largest 
asteroid, although low significance prograde rotation persists even 
for smaller asteroids. 



4.2 Spin- up of protoplanets 

The numerical experiments presented in |3] show that ac- 
cretion of particle sizes of a few cm to a few m and a 
solids-to-gas ratio around unity or higher (see Figure 3) 
leads to a systematic prograde rotation of protoplanets, 
in agreement with the observed trend in the largest as- 
teroids. This moderately high solids-to-gas ratio points to 
the sedimentary mid-plane, or to den se particle clumps ap- 
pearing in such a mid-plane layer (jjohansen et al.l 120071 ; 
ICuzzi. Hogan. fc Sha rif! 2008), as a prime site for the growth 
of protoplanets. The evaporation of gas by irradiation from 
the central star or neighbouring O stars leads to a simi- 
lar i ncrease in the influenc e of the dust particles on the 
gas jThroop fc Ballvl 120051; lAlexander fc Armitagel 120071 ; 
iJohansen. Youdin. fc Mac LowH2009r i. 

Our simulations predict that the largest Kuiper belt 
objects should preferentially rotate in the prograde direc- 
tion around their axis. A few exceptions to this trend are 
objects that have suffered massive impacts (e.g. Haumea, 
iBrown et ai1l2007r i. 

Pure accretion from a circumplanetary disc would lead 
to a protoplanet rotation period of only a few hours, near the 
centrifugal break up speed. However, accretion from a cir- 
cumplanetary disc need not be the only mass source. Growth 
phases with little or no accreted angular momentum can oc- 
cur when the protoplanet accretes isolated dust particles 
and pebbles or km-sized planetesimals or other protoplan- 
ets. Protoplanets thus need to obtain between 10% and 50% 
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Figure 9. Cumulative spin period distribution of Kuiper belt 
objects and asteroids larger than 200 km in radius. Most objects 
lie in a narrow 5 < P < 20 hour period range irrespective of their 
other physical properties. Outliers Haumea (H) and Pluto (P) are 
marked; they are believed to have suffered massive collisions. 



of their mass from the accretion of macroscopic dust par- 
ticles in a gaseous environment. The lower limit is relevant 
if all pebbles are accreted at the full Keplerian angular mo- 
mentum, while the upper limit is relevant if 80% of the mass 
accretion occurs through pebble sandblasting or through ac- 
cretion of "planetesimals" with high T{. 

The sub-Keplerian sandblasting effect described in £|3 . 2 1 
means that the mechanism for obtaining prograde rotation 
works best in particle-dominated flows or in regions of low 
radial drift. The radial drift can be reduced or even re- 
versed in persistent pressure bumps appearing in the gas 
l|Haghighipour fc Boss! 120031 ) . Recent studies have shown 
how planetesimal form ation by self-gravity dJohansen et al. 
20071 ; iKato et all l2009|j and by coagulatio n (|Kretke fc Lin 
20071 ; iBrauer, Henning. fc Dullemondl fcoOS) proceed in such 
convergence regions in the flow of solid particles. 



4.3 Spin-up in absence of gas? 

It was recently proposed by [Schlichtin g~fc Saril (|2007l j that 
inelastic particle collisions will lead to accretion of prograde 
rotation from a circumplanetary particle disc during the gas- 
free stages of planetary growth. Such a mechanism may be 
able to provide additional prograde angular momentum af- 
ter the giant impact stage. This is especially relevant for the 
terrestrial planets which are believed to have been accumu- 
lated over 10-100 m illion years in a gas free environment 
ijLunine et al.ll2010T i. Mercury, Earth and Mars all have pro- 
grade spins, but the spins of the terrestrial planets have been 
substant ially modified by long-te rm perturbations by other 
planets (ILask ar fc Robutell 19931 ). 

Focusing instead on the largest asteroids and Kuiper 
belt objects, our numerical simulations show that pebbles 
and rocks are accreted very efficiently by the protoplanets. 
In the main numerical simulation presented in [J3]the proto- 
planet+disc system accretes 1/7 of the protoplanet mass in 
five orbits (~ 25 years at 3 AU). Such high accretion rates 
are facilitated by the presence of gas in the Hill sphere. This 
indicates that young protoplanets can accrete lots of mate- 
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rial already during the protoplanetary gas disc phase. The 
actual ratio of accretion during and after the gaseous disc 
phase is nevertheless relatively unconstrained and deserves 
further attention in the future. 



5 CONCLUSIONS 

We present computer simulations of the accretion of peb- 
bles and rocks onto young protoplanets in a gaseous nebula 
environment. Our findings can be summarised as follows: 

• Particle streams interacting with the gas far out in the 
Hill sphere of the growing protoplanet dissipate enough en- 
ergy that the particles spiral down towards the central object 
and form a prograde accretion disc. The prograde motion is 
transferred to the protoplanet as the particles are accreted. 

• In order to explain the observed rotations of the largest 
asteroids and Kuiper belt objects, we propose that proto- 
planets accreted a major fraction (10%-50%) of their mass 
from pebbles and rocks during the gaseous solar nebula 
phase. 

• A specific prediction of our model is that most of the 
largest Kuiper belt objects should possess prograde spins. 

Since the Keplerian surface frequency depends only on 
the density of the protoplanet, accretion from a Keplerian 
particle disc explains how asteroids and Kuiper belt objects 
can have similar rotation periods even though the orbital 
time-scale is much longer in the Kuiper belt. Our results 
imply not just a universality in the rotation period of solid 
bodies in the solar system, but also in the formation process 
and in the important role of gas and pebbles for the growth 
of hundred kilometre scale protoplanets. 

Future numerical work should focus on the extension of 
the presented model to include turbulence driven by global 
instabilities and more realistic solid physics, such as a dis- 
tribution of particle sizes, collisional viscosity and collisional 
fragmentation. 
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APPENDIX A: DYNAMICAL EQUATIONS 

We simulate the couple d motion of particles and g as in 
the shearing sheet frame (|Goldreich fc Tremainelll980l ). The 
frame rotates at the Keplerian frequency J?k at an arbitrary 
distance ro from the central star. Axes are oriented such 
that x points radially away from the star, y points along the 
rotation direction of the disc, while z points vertically out 
of the orbital plane. The particle component is treated as 
Lagrangian particles, each with a position vector x = (x,y) 
and a velocity vector v = (v x , v y ), with dynamical equations 



x — 



v 



x , 



y — v y - (3/2)i7 K x , 
v x — +2{2kv v — 
1 



_ GM p x _ JL_ 
(x 2 + y 2 ) 3 ' 2 ~ n {Vx 



_ GM p y 

2 UKVX (x 2 +y 2 ) 3 / 2 



Tf 



-(Vy ~ Uy) . 
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(A2) 
(A3) 

(A4) 



We measure the velocity in the y-direction relative to the 
linearised Keplerian shear u y ^ = —(3/2)Ok%, which gives 
rise to a change in the Coriolis force appearing in equa- 
tion ()A4[) . An additional force contribution from the pro- 
toplanet of mass M p fixed in the centre of the coordinate 
frame is added. We also consider the drag force from the 
gas, acting on the friction time-scale Tf and proportional to 
the velocity difference between particles and gas. The gas 
velocity u = (u x ,u y ) at the position of a particle is found 
by spline inter polation from the nearest t hree grid points in 
each direction (|Youdin fc Johansen! [20071 ). 

The evolution of gas is solved on a fixed Cartesian grid. 
The equations of motion for the gas velocity field u, relative 
to the Keplerian shear, and the continuity equation for gas 
density p g , read 
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(A6) 
(A7) 



In equation (|A5|) we include the possibility of adding a radial 
pressure gradient d In P/d In r through the sub- Keplerian ve- 
locity 

Av = -~ — — c B . A8 

2 r dmr 

We use An = except i n H3.2I where we set Av = 0.05c s 
|Nakagawa. Sekiva. fc Havash i 1986). The gas also feels the 
gravity of the protoplanet. For the back-reaction friction 
force from the particles we calculate the local solids-to-gas 
ratio through e = p p /Pg- 

The dynamical equations for the evolution of particles 
and gas appear in 3-D as 
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Here we subject particles to the vertical gravity component 
of the central star, but ignore the same term for the pressure- 
supported gas. 



APPENDIX B: CODE VALIDATION 

To validate the implementation of the protoplanet gravity in 
the shearing sheet frame, we test here t he accretion of test 
parti cles against other published results |Schlichting fc SarH 
12007? ). 

We consider a protoplanet of dimensionless mass 
GMp = W^n^H 3 in a 2-D box of size L x = L y = 0MH. 
Here Ok is the Keplerian orbital frequency, while H is 
the gas scale height. The Hill radius of the protoplanet is 
-Rh ~ 0.0077?. We place 10,000 particles randomly in the 
box, initially excluding particles from the region within the 
Hill sphere. The particle positions and velocities evolve in 
response to the Coriolis force, the tidal force from the cen- 
tral star and the gravity of the embedded protoplanet. For 
this test we ignore drag forces between particles and gas. 

As the particle positions and velocities evolve, we track 
the first crossing of various distances from the protoplanet 
and calculate the mean specific angular momentum per par- 
ticle for all the accretion shells. The results are shown in 
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Figure Al. Code validation test. The left plot shows the mean specific angular momentum, (£z), of test particles passing circles of 
various radius -Race (x-axis). In the limit i? aC c 3> i?mn the slow prograde rotation of the nebula is transferred to the "planet" [the dotted 
line shows the analytical expectation £ z = (l/4)i?K-R 2 cc ]j but particles rapidly decrease their angular momentum as they travel into the 
Hill sphere. The dashed line indicates prograde Keplerian specific angular momentum. The right plot shows the ratio of accreted angular 
momentum to the prograde Keplerian value. These non-interacting planctcsimals induc e a net retrograde rota tion on the protoplanet. 
The measured angular momentum flux is in good agreement with results published in ISchlichting fc Sari ( 2007]) ■ 



Fig. IA1I For accretion outside of the Hill sphere, the motion 
is prograde, fo llowing the trend l z = (1/ 4) Qn R^rr derived 
analytically in iDones fc Tremaind j 19931 ). However, deeper 
in the Hill sphere the net angular momentum is approxi- 
mately zero. Particles transported to the vicinity of the pro- 
toplanet by the Keplerian shear lose their angular momen- 
tum on the way and induce slow retrograde rotation to the 
protoplanet. The nu merical result shown in F ig. IA1I is in 
good agreement with lSchlichting fc Saril l|2007l ). 



APPENDIX C: PARTICLE SIZES 

We characterize particles in terms of their dimensionless 
friction time Qkti (where J?k is the Keplerian orbital fre- 
quency). The translation to a particle size depends on the 
assumed disc model. Small particles are subject to Epstein 
drag, valid when the particle radius is smaller than the mean 
free path of the gas (see below). The friction time in the Ep- 
stein regime is given by 



(E P ) _ ftnp»a 



p. a 



PgC s 



(CI) 



where the second step applies in the disc mid-plane. Here 
p. is the material density of the solids, a is the radius of 
a solid body, p g is the gas density, c s is the sound speed, 
while S g is the column density of gas. Epstein drag depends 
on gas density, but in practice the gas density fluctuations 
are negligible (order 1%) in our subsonic flows, so we ignore 
them. Solving for particle size gives 
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Applying a simulation with given n values to different 
locations in the disc with different gas column densi- 
ties changes the physical particle sizes corresponding to a 
given value of Pu ff. We follow the Minimum Mass So- 
lar Nebula model dWeidenschillin j Il977bl ; lHavashil Il98ll ; 
ICuzzi. Dobrovolskis. fc Champnevl 19931) . with column den- 
sity X" g and gas scale height H given by 



H(r) 



1700 gem" 



/ r \ 5 / 4 
0.033 AU (-) 
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(C3) 
(C4) 



Considering, as in the main text, r = 3AU and E s = 
300gcm~ 2 or r = 10 AU and E g — 50 gem -2 , yields particle 
sizes of a = 6 cm and a = 1 cm, respectively, for HkT{ = 0.1. 

The Epstein regime of free molecular flow ceases to ap- 
ply once the particle radius exceeds (9/4 of) the gas mean 
free path, 
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where fi = 3.9 x 10~ 24 g is the mean molecular weight and 
(j mo i = 2 x 10~ 15 cm 2 is the molecular cross section of molec- 
ular hydrogen l|Nakagawa et al.1l 19861 ). Epstein drag applies 
as long as 



2 £ g <T mo \ 



300 gem" 



2 g cm" 
H/r 



0.043 / V3AU 



(sau) ■ ^ 



Thus Epstein drag is the relevant regime for the application 
of our model to particles with Qktj < 1 beyond r = 3 AU. 
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Figure El. Convergence in grid resolution. The plots show the accreted angular momentum for runs with 256 2 grid points (left panel), 
512 2 grid points (middle panel) and 1024 2 grid points (right panel). In all cases a net positive angular momentum is accreted, although 
periods of retrograde accretion occur as well. The inner radius of the simulations, where particles are considered to be accreted to the 
protoplanet, is reduced proportional to the resolution. 



APPENDIX D: PARTICLE COLLISIONS 

In this section we argue that particle collisions are negligible 
compared to gas drag in most of the simulation domain. The 
time-scale for collisions between particles is 



(Dl) 



where c p is the local velocity dispersion ("temperature") of 
the particles and A p is the mean free path for particle col- 
lisions. The mean free path in a granular medium can be 
calculated from the particle number density n p and colli- 
sional cross section <j c , 



1 
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Assuming spherical particles of radius a, this yields a colli- 
sion time-scale 



ap. 



(D3) 



3c p p p 

The ratio of Epstein friction time to collision time is (ignor- 
ing factors of order unity) 

21 ^ c p p p 

T c C s f}g 
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Particle collisions become important when p P /p s > c a /c p , 
i.e. for flows with large particle loading. For the typical 
value c p ~ 0.01c s , particle collisions are important when 
p P /pg > 100. Thus we can ignore particle collisions in the 
average state of the simulations presented in the main pa- 
per. Nevertheless particle collisions will become important 
in the dense regions close to the protoplanet, and indeed 
particle collisions may be important for efficient accretion 
of circumplanetary material onto the protoplanet. 



have run lower resolution (256 2 grid points, 250,000 par- 
ticles) and higher resolution (1024 2 grid points, 8,000,000 
particles) variations of the main simulation. We set the in- 
ner accretion radius to 4 grid points in all cases. The ac- 
creted angular momentum is shown in Fig. IE1I It is clear 
that the behaviour at all three resolutions is highly time 
variable, with periods of prograde and retrograde accretion 
interspersed. However, in all cases a net positive angular mo- 
mentum is accreted, inducing between 25% and 50% of the 
surface Keplerian frequency to the protoplanet. There is a 
trend for accreting stronger prograde rotation with increas- 
ing resolution. 



APPENDIX F: DEPENDENCE ON BOX SIZE 

We proceed to test the dependence of our results on the size 
of the simulation domain. For this purpose we ran a low 
resolution simulation with the radial and azimuthal extent 
expanded by a factor 2.5, to L x = L y — 0.1. We kept the 
grid and particle resolution constant by using 640 x 640 grid 
cells and 1,562,500 particles, corresponding in resolution to 
the 256 2 simulation shown in the left panel of Fig. IE1I The 
accreted angular momentum is shown in Fig. IF1I The larger 
box size has a net positive effect on the prograde rotation 
of the protoplanet. The sign of the angular momentum is 
more consistently positive than in the small box, and there 
are fewer periods of retrogade accretion. 

This paper has been typeset from a TpX/ ETpX file prepared 
by the author. 



APPENDIX E: CONVERGENCE TEST 



To quantify the effect of grid resolution and inner accre- 
tion radius on the results presented in the main paper, we 
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Figure Fl. Convergence in box size. The plots show the accreted 
angular momentum for low resolution runs with L x = L y = 
0.04H (top panel) and L x = L y = Q.IH (bottom panel). Pro- 
grade motion is more readily accreted to the protoplanet in the 
larger box. 



